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ABSTRACT 


A new wind-tunnel technique for measuring various aerodynamic 
derivatives of an aeroelastic model is presented. The technique applies 
free- flight procedures to a model flown in the wind tunnel on the two- 
cable-mount system. The complete equations of motion are presented. 

In the case of longitudinal motion, it is theoretically possible 
to uniquely determine each of the aerodynamic derivatives by measuring 
the model response to a steady-state sinusoidal oscillation of the 
horizontal tail. A comparison between free- flight and wind-tunnel equa- 
tions shows that, due to the added mount system restraints, the equa- 
tions can be solved for each derivative uniquely. However, introduction 
of an error into the model response investigated caused the solution to 
become ill-conditioned, resulting in equations similar to those used for 
determining the derivatives in free flight. 

In the lateral equations of motion, a basic free- flight assump- 
tion of single-degree- of- freedom response in roll allows the experimental 
verification of the dynamic approach to derivative measurements. Experi- 
mental results obtained on an aeroelast ically scaled model, tested both 
statically and dynamically in the wind tunnel, verify the application 
of this new testing procedure. 
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CHAPTER I 


INTRODUCTION 

In the study of aircraft stability and control a knowledge of the 
aerodynamic characteristics of the aircraft is of basic importance. 

The aerodynamic characteristics of a flexible airplane are difficult to 
measure or predict analytically, especially at transonic speeds. The 
purpose of this thesis is to present a new wind-tunnel testing technique 
for the measurement of certain aircraft aerodynamic characteristics 
known as the aerodynamic stability derivatives. An aerodynamic 
stability derivative indicates the rate of change of a force or moment 
acting on an airplane with the motion or variable causing the force or 
moment . 

Structural components of an aircraft are often flexible enough to 
be considered as nonrlgid. That is, such aeroelastic phenomena as 
flutter, divergence, and control surface effectiveness must be investi- 
gated. It is also realized that structural flexibility may have 
appreciable Influence on the aerodynamic derivatives and thereby affect 
the overall flying qualities of the aircraft. Aerodynamic considerations 
of structural flexibility occurring in the transonic speed region can be 
especially difficult since no dependable aerodynamic theories are 
available. Consequently, the necessity for measuring the aerodynamic 
derivatives of a flexible aircraft is apparent. 

Within recent years a mount system has been developed (Ref. 4) 
which permits the "free-flight" behavior of an aircraft to be simulated 
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in the wind tunnel. This mount system, referred to as the two-cable 
mount, is shovn schematically in Figure 1 (p. 7 ). The model is held 
by two mutually perpendicular cables passing through pulleys in the 
fuselage and attached to the tunnel walls. The cables are kept under 
tension by stretching a soft spring in the rear cable. Remotely 
operated trim controls are Ja^a^Lded on the model. The mount system was 
originally designed for the testing of aeroelastic effects such as 
flutter. Since the equations governing motion on the mount system are 
quite similar to those of flight, it was soon realized that the mount 
offers a potential for measuring aerodynamic derivatives of a flexible 
model in the wind tunnel. 

The equations governing model behavior in the wind tunnel are 
essentially the free-flight equations modified by the addition of mount 
system restraints. It is therefore possible to apply test techniques 
similar to those presently used to obtain free-flight data. The tech- 
nique selected for investigation Involves measuring the model response 
to a sinusoidal steady-state excitation provided by the model control 
surfaces . 

The equations governing model motion on the two-cable-mount system 
are derived in terms of model mass properties, linearized mount 
restraints, and unknown aerodynamic derivatives. Since fore and aft 
motion is not provided by the mount configuration analyzed, the 
equations of motion are presented In five degrees of freedom. The 
equations are simplified by separating them into two independent groups. 


V 
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The longitudinal equations include vertical translation and pitching 
motions. The lateral equations include side translation, roll, and 
yawing motions. Each of these sets of equations is treated separately. 

A technique for measuring the longitudinal aerodynamic deriva- 
tives is presented. The longitudinal equations are expressed in two 
degrees of freedom. A sinusoidal steady-state forcing function generated 
by the horizontal tail is introduced into the equations of motion. A 
steady-state sinusoidal response is assumed, and the resultant equations 
are expanded in terms of measurable model response and the unknown 
stability derivatives. It is assumed that the derivatives are indepen- 
dent of frequency; therefore, the equations of motion are valid at each 
discrete excitation frequency. From the model response measurements at 
several frequencies, a set of redundant equations is generated which 
can be solved for the unknown derivatives. A least-squares solution is 
used to obtain the derivatives from the set of redundant equations. 

Since no experimental data are available for longitudinal motion, 
a numerical example is given to determine the effect of measurement 
errors on the derivatives. Equations for the two- cable mount show that 
it is theoretically possible to separate the results into uniquely 
determined stability derivatives. Greenburg (Ref. 2) states that results 
obtained by the dynamic technique from free- flight measurements appear as 
linear combinations of the aerodynamic derivatives caused by a dependency 
on the aircraft response. A comparison between flight and wind-tunnel 
equations shows that the mount system restraints allow separation of 
the unknown derivatives. In practice, however, the magnitude of thiB 
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restraint term is overshadowed by small errors in response measurements, 
which result in an ill-conditioned problem when solving for each of the 
derivatives uniquely. 

The problem of determining each of the aerodynamic derivatives 
appearing in the lateral equations is more complex than that of the 
longitudinal case because of the added degree of freedom and its asso- 
ciated derivatives. However, Etkin (Ref. 1) states that for many 
conventional airplane configurations the roll equation in flight can be 
simplified and treated as a single degree of freedom. 

In order to verify the basic dynamic approach to derivative 
measurements, an experimental technique for measuring aileron effective- 
ness and damping -in-roll stability derivatives is presented. Modified 
flight techniques are applied to the wind-tunnel tests . The experimental 
technique Involves the measurement of model response to a sinusoidal 
steady-state forcing function generated by the ailerons . In order to 
satisfy free-flight requirements (single degree-of -freedom response), a 
parametric study of the two-cable-mount system was run in an effort to 
force the model to behave, essentially, as a single degree of freedom 
in roll. If the aerodynamic derivatives are independent of frequency, 
the roll equation is valid at each discrete aileron frequency. Measuring 
the model response as a function of aileron frequency generates a set of 
redundant equations which are solved for the unknown derivatives. In 
this manner both the aileron effectiveness (Cj^) and damping-in-roll 
derivative (C^) are determined. 
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Experimental results are presented for a l/l9-size aeroelastic 
scaled model of a large, subsonic, multijet. cargo airplane. Experimental 
results are obtained dynamically, as described, to determine and 

C, and statically to determine C 7 . A brief description of the static 

4 p 4 6 

mount and test procedure is given. A comparison betveen static and 
dynamic tests shows good agreement within the basic assumptions made. 

Based on these results, the application of flight techniques to 
scaled models flown in the vlnd tunnel on the two-cable mount offers a 
potential for making quantitative measurements of the effect of flexi- 
bility on aircraft stability derivatives. This new testing procedure 
offers the aircraft designer early estimates of the stability and control 
characteristics of future aircraft configurations. 



CHAPTER II 


EQUATIONS OF MOTION 
2 . 1 Introduct ion 

In this Chapter the equations of motion governing model behavior 
on the two-cable-mount system are developed. Cable restraints are 
presented as stiffness influence coefficients. The assumption of small 
perturbations from trimmed flight makes the equations of motion linear 
and allows them to be separated into longitudinal and lateral degrees of 
freedom. 

The mount configuration shown schematically in Figure 1 is 
analyzed. The model is held by two cable loops ; the upstream cable is 
in the vertical plane, and the downstream cable in the horizontal plane . 
A soft spring in the rear cable keeps the system under tension. The x, 
y, z axes form a fixed right-hand coordinate system with its origin at 
the model center of gravity. The x axis is directed upstream in the 
tunnel, and the z axis is in the direction of gravity. The equations 
of motion are limited to five degrees of rigid-body freedom since fore 
and aft motion is not provided by the mount configuration analyzed. 

Assuming positive displacements as shown in Figure l, the dynamic 
equations of motion can be written as follows: 

Vertical translation: 

Summation of forces in the z direction 

Z A + Z c + mg - mi! (la) 
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pitch: 

Summation of moments about the y axis 
M A + Me « Iy8 

Side translation: 

Summation of forces in the y direction 
Y A + Y C - my 

Roll: 

Summation of moments about the x axis 
L k + h m hP- J XZ* 

Yaw: 

Summation of moments about the z axis 
N A + N C " h* ' X XZ# 

where 


(lb) 


(lc) 


(Id) 


(lc) 


z * vertical translation of model center of gravity 
0 ® rotation about y axi6 

y ■ lateral translation of model center of gravity 
0 ■ rotation about x axis 
if > rotation about z axis 
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and 

L,M,N are moments about the x, y, z axes, respectively 
Z,Y are forces along the z and y axes, respectively 

The subscripts A refer to aerodynamic forces and moments; the 
subscripts C, to forces and moments generated by the mount restraints. 

Etkin (Ref. l) shows through the use of small perturbation theory 
and other assumptions that the complete set of eolations of motion In 
flight can be separated Into two Independent groups. These tvo groups 
are referred to as the longitudinal and lateral equations of motion. 

The equations for longitudinal motion Include fore and aft motion, 
vertical translation, and pitch. The equations for lateral motion 
Include side translation, roll, and yaw. The basic longitudinal, and 
lateral equations of motion for the two-cable-mount system will be 
presented in the following sections . 

2.2 Derivation of Equations of Motion 
2.2.1 Aerodynamic Force Components 

2.2.1 .1 Longitudinal Forces .- Lift and drag forces (L and D) 
are defined to have directions normal and parallel, respectively, to the 
relative wind vector V, as shown in Figure 2. We have the following: 

V - % + z 
• 

a ■ 8 + - 
U 

o- te 

e - e t + e 


for z < < U 
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where 

V = relative wind vector of the airplane 
V 0 = wind vector along x axis 
U = wind-tunnel flow velocity 

- pitch angle for trimmed flight 
0 = pitch angle perturbation about trim 

Vortical translation: 

Referring to Figure 2, the sum of the aerodynamic forces in 
the z direction results in 

Z A = -L cos (a - ?) - D sin(a - 0) 

Since it is assumed that z < < U, then 

Z A = -L - D(a - 9) (2) 

The aerodynamic forces are normally expressed in terms of their coeffi- 
cients 


L - qSC L 
D » qSC D 


The nondimens ional coefficients, and C D' represent the lift and 

drag generated on an airplane for a given dynamic pressure q and 
representative wing area S. ^Dynamic pressure q is defined as 
q * ^ pU 2 , where p is the testing medium density.! 


¥ 
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Equation (2) can be vrltten as 


Z A = -qs 


+ C D u] 


Making the basic assumption that the lift coefficient is a 

linear function of angle of attack a and tail angle 5, ve can 
express as a Taylor series about the trim point 8^ . therefore, 


oCr 3Cr 

C T = C T + — * Ax + — ^ AS 

L ^ 3a 38 


vhere 

= lift coefficient at trim 8^. 

Aa = a (angle of attack perturbation) 
AS = 6 (tail angle perturbation)'*' 


3Ct 


Therefore, equation (3) can be vrltten as 


“ Z Aq " qS + C D uj " qSC L6 5 


‘‘'Sign convention used is +8 is a trailing edge up vhich 
generates a negative lift. 


co|ro 
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Z A - \ 


**[ c J e+ §) + c »j] 


Zp^ - trim aerodynamic forces 
Z A = -qSCj^S 

Pitching motion: 

The aerodynamic pitching moment is normally expressed as 


% = <jScCn 


C m is a nondimens ional coefficient representing the moment generated 
for a given dynamic pressure, wing area, and a representative length c. 
(c is the mean aerodynamic chord. 1 ) Seckel (Ref. 5) shows that the 
aerodynamic moment can be expressed in a Taylor series expansion about 


0 t as 


dM A dM A . dM A • dM A 

M. = M. + + 

A *o da ad ae as 


1 The mean aerodynamic chord is defined by the formula 

r b/2 

/ C 2 dy where C is the local chord, b is the span, S is 


the wind area, and y is the lateral coordinate. 
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where 

s static trim condition 
Aa = a, angle-of -attack perturbation 
Ax - a, rate of angle-of -attack perturbation 
A9 = 8, rate of pitch-angle perturbation 
A5 =8, tail-angle perturbation 

Substituting equation (5a) into (5b) results in. 


M a = M. 


qSc - — a + — — a + — 7 — 9 + — — 8 I 
da da d0 d6 J 


( 6 ) 


In order to handle the aerodynamic terms more conveniently, it is 
common practice to nondimensionalize the aerodynamic derivatives in the 
following manner: 



Substituting the preceding terms into equation (6) and expanding the 
expressions results in 
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ma*%*^[v + SxMjSlMj (c w + v 4 ] + ^ 


( 7 ) 


where 


M A = < i SSC m 6 b 


(7a) 


2.2 .1.2 lateral Forces .- Hie lateral aerodynamics are determined 
in a manner similar to that in which the longitudinal aerodynamics are 
derived. The lateral aerodynamic equations can be written as follows 
(Ref. 4): 

I* A = ^SbCj 
N a = qSbC n 
Y A = qSC y 

where 


Cj ■ rolling -moment coefficient 
C n = yawing -moment coefficient 
Cy = side-force coefficient 

lhe representative length for the lateral equations is the wing span b. 
Hie angle of sideslip 3 is defined as ^ - ijr. Hie lateral aerodynamic 
equations can now be written as follows: 


L A - 1Sb(c, B § + |g i - c, ( t + C, r lo ♦) * L A 


( 8 ) 
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"a • ««*•(% g ♦ \ k* - %*♦ \ is*) + *A 
«*(<% - C D>C + <X« I + C y p 3} i - V + °y r W ♦] + f A 


(9) 


(10) 


where 


l a = 



*A = qSbCy^Aj 


6 A - aileron angle perturbation 

The expressions for L A , N A , and Y A are simplified forms corresponding 
to a fixed rudder configuration which is utilized in wind-tunnel testing 
procedure. The aerodynamic derivatives are nondimenslonalized as 
follows : 


dC 7 

se ‘e 

dOj b 

W = 

‘ " 0 t_ 

| 2U *p 

|£l R 4£ir - 
aF su 7F\ 

i- C 7 
2U *r 


\ 2U) 


S \ 2 U/ 


. 

d g n m b_ dCp 

b „ 

^n b ^n 

b 

ST’S 

" 2 U *(&) 
\zu) 

“ * Cn- 

| 2U °P 

41 '"SHT 

2U Cn *“ 




M 



17 



d°y _ b_ dC y b 
d$ " 2U ^jb j " 2U 



df 2U Jjb 
V2U 



dC 7 

^ ■ s 




dCy 

* 6 A 



2.2.2 Mount Restraint Force Components 

The two -cable-mount configuration used during this investigation 
is presented in Figure 3» The linearized cable restraint forces for the 

longitudinal degrees of freedom are derived in Appendix A. Equations 
are developed to determine the longitudinal spring constants in terms 

of mount system geometry, tension, etc. Hie longitudinal and lateral 
cable restraint forces can be expressed as follows: 

z c * " K ze 0 " K zz z (u*) 


“ *c 0 ” ~ 

- -Kj^yy - Kflfj - 

Me * “fyyy “ 



Y C • -Kyyy - - Ky t *r 


(lie) 


nit w »Bi ftm f 
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where K-y Is the spring constant relating a force or moment In mode i 
due to a displacement in mode J . 

For the configuration under investigation (Appendix A) the stiff- 
ness influence coefficients are as follows: 

Longitudinal stiffness Influence coefficients 

+ 5fJ 

"ze ' 2 [ a ^ ‘ £ h “ e h 8ln h ' 5* e c «%] 


z c 0 " " K ze e t 
Kqz = K ze 


Tr 


he * 2 ‘‘ if * 2 ‘ t r «• h * ^ sin Pp + e cos Pp 
+ ^(e cos Pp + h sin pp)*j 


Lateral stiffness influence coefficients 


% ’ 2 % * 2 1 ; cw2e » 


Ktf • 0 
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k Y* “ 2 e - 2 cos 3 _(d sin P R + a cos 0 R ) 
v ^ Lr ‘ 

^Y * " 0 

*# - + ,in 

Y ~yf> = = 0 

k *y = 

- 2 cT f(§: + 008 Pf) + ^r ] 8 cos h + d 8in Pr 
i 2 I 

+ ^(a cos Pr + d sin 0 R ) J 

2.2.3 Complete Equations of Motion 

Having developed the aerodynamic and mount forces, equations ( 4 ), 
(?)> (8), (9) (10), (11a) through (lie) are substituted into 

equations (la) through (le) resulting in the following equations of 



motion: 


Vertical translation: 


Z. + 2. 

A Aq 


,S (%( 9 * u) + C D u] * Z C 0 • V - “zz 2 * " 8 ■ ■* 


(12a) 


Pitch: 


M A + % + 9 s * (v + C %S + lj C HlU + Ij (C n>i + V®] 


+ «Co ” Kez 2 “ ^e 6 ■ x v® 


(12b) 


Side translation: 


\ * ’ s (<% - C D>5 V * Cy r A *] 


v - v - V 


my 


(12c) 


Roll: 


L A + 


qSb (% U + c lp Zjt ’ <V + c h ZB ♦] ' *fr y ’ 


The longitudinal equations can be simplified further by considering the 


static trim condition 0 t . At 0 t let 0»0*z«z-z-O; 8»6 0 . 

Therefore, 


-qSC L& 5 Q + Z Aq + ^ + mg - 0 
qScCmgSo + MAo + MCq * 0 


Static trim equations 


The dynamic longitudinal equations of motion for small perturbations 
about the trim point 0^ become the following: 

Vertical translation: 




Kze® ” K ZZ z * nz 


Pitch: 


(13a) 


M a + qSc 


(v + c % 5 * #j Si I * 3J ( Sx * W 9 ] • • “ee 9 ‘ ^ 
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The equations of motion are linear second-order differential 
equations with constant coefficients. The analysis does not include 
the dynamic properties of the mount itself such as cable inertia and 
pulley damping. A simplified approach is to neglect cable inertia and 
add viscous damping terms proportional to displacement rates in each 
mode. The terms C^i, C 0 0, Cyy > C^\|r, and C^J are added to each of 

the equations of motion, respectively. Where 

PS 

z = V“m~ 

= 

“W = 

°%r 

The damping factors £ z , and ^ are estimated for the 

pulley configuration being investigated. Substituting these terms into 
equations (12c), (12d), (12e), ( 13 a), ( 13 b), and rearranging, results in 




C Z = 2 ^Z ma ^Z 

c e = 2£ei Y <jfce 

Cy = 2^ Y maA fY 

S = 2 ^ I z a H 


the following equations: 


2k 


Vertical translation: 


- * [S ^ * '!>> + 2Cz mCt zzj^ + *2z z + + tfCjJe • 


Pitch: 


f® " ftr ( Si + <V • 2 Vy %)] 4 + (Kee " <^V e 


qSc^ _ .. qSc 


- -3- Si 2 - V +K ez z - % (l^b) 


Side translation: 


C D ) - 25 yfflaY 3 fJy + Kyjy - |jj- Cyjt + (Ky0 - qSC L )0 


■ lr V + (qS( Vs + % )+ = y a 


(14c) 


^ - (lr c . p - ‘^-v-r v 

* (?SbC, ♦ Kjfr)* - ^ C, y + K^y * L A 


(ltd) 
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Yaw: 

V - C n r - 2 ?*iz>v)i ♦ («»% ♦ V* - ^ S> 

+ *«» - W - ^ ♦ V * *A 

(l4e) 

Equations (l4a) through (l4e) completely describe, within the 
framework of the assumptions made, both the longitudinal and lateral 
motion of the model in the wind tunnel. The equations are expressed 
in terms of model mass properties, tunnel test conditions, mount 
properties, and aerodynamic derivatives. 


CHAPTER III 


TECHNIQUE FOR MEASURING LONGITUDINAL 
AERODYNAMIC DERIVATIVES 

The equations of motion governing model behavior in the wind 
tunnel are quite similar to the free-flight equations modified by the 
addition of mount system constraints. Hence, it is possible to apply 
test techniques similar to those used to obtain aerodynamic derivatives 
from free-flight tests. The technique selected involves measuring the 
dynamic response of the model to a sinusoidal excitation. In this 
thesis, the excitation is provided by a sinusoidal movement of the 
horizontal tail. Greenburg (Ref. 2) presents this technique for the 
free-flight case. 


3.1 Derivation 

For simplicity, it is assumed the pulley damping and cross - 
coupled mount stiffnesses in equations (l4a) and (l^b) are zero (i.e., 

Cg * ^ = Kqz = K 2 Q =0). The equations of motion become 

mz + y (Cj^ t Cjj)z ^ Y uyrj z + qSCy^Q ~ (15a) 


V - + <*68 ' 


Sr v 


* ■ Sa (15b) 
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In order to provide a dynamic forcing function, let the tail 
deflection be 6 = 6 0 e icut . A sinusoidal deflection of the tail will 
provide translation and pitch forcing functions of the form: 

Z A (t)« -qSC^Boe 10 * (l6a) 

M A (t) = qgcC^Boe 1 ^ (l6b) 

Substitution of the forcing functions in the form Z A (t), M A (t) 
into equations (15a) and (15b) results in. 


mz 


♦ f (V ♦ 


Cjj)z Kyy z + qSCr Q » 


A® 


iajt 


(17a) 


2 __2 
M * (“09 - 9S5<^)9 - y- C^i 


‘ IF V = 9Sc<V> 0 ei«* 


(17b) 


For sinusoidal motion, the steady-state response becomes: 


z . z 0 e 1( “* 4(l l > 

i = 

z = z o e( aA+ 0l) 


8 = 

8 » loeoe 1 * 0 *^* 

0 = -o^0 o e i ( a5t+ ^2) 
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0^ and 02 the phase angles between tall and model displacements. 
Substituting the response functions Into ^ITa) and (1Tb) results In the 
following equations: 

Vertical translation: 

+ Cj^+qS&Q) = (mo? - Kgzhoe^ 1 ( l8 *) 

Pitch: 

S [ '4 «**] ♦ **[- ^ ^ 

+ J^qS50 o ei ^ - iai ^ z^ 1 ^ + c m 8 G' qS ^ 5 q3 

= (I yO? - Kqq )*^ 2 

(18b) 

If the aerodynamic derivatives Cj^, C D , Cj^, C^, C^, 

and are assumed to be independent of frequency, each of the equa- 

tions of motion can be used to generate a set of redundant equations 
which can be solved, using a least-squares method of solution (refer to 
Appendix B) for each of the aerodynamic derivatives. The least -squares 
method of solution presented in Appendix B is a mathematical procedure 


ij 


* 
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for solving a set of redundant equations. Equation (l8a) can be expressed 
as a function of frequency in the following form: 


t qSQQe 1 ^ 2 +ia3y ZqC 1 ^ 1 




J ~ 1,2,3, . . . N 

(19a) 

0 O , z o> $l> and $2 31,6 measure<i from ®odel response as a function of 

the tail frequency go. Similarly, equation (l8b) can be written as 



= Qiya£ - Kee) 0 o ei ^ 2 ] 

0J=0)J 


J " 1*2,3* • • • N 


(19b) 
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Equation for vertical translation: 

Equation (19a) can be rewritten in the form, 

+ A J2 c D + A J5 C Ie = b J ( 19c ) 

where 

A J1 - + io) ^ e^lj 

J C0=C0j 

A J2 - * f L^ 1 ] 

CJ=COj 

A J3 = + < 3S8 0 

bj = [W - K Z z)z 0 e 1 ^ 1 ] J = 1,2,3, • • • N 

O^CDj 


Applying a least-squares solution to equation (19c) results 
in the following (Appendix B) : 


*IM 2 

N 

Z^l^+AjlAJg) 1 

N 

Z (A Jl A J3 + A Jl A 53 ) 


J=i 

J-l 

J-l 


N 

Z (A 5l A J2 

N 

+ A J1 A J2) 

N 

C D 

J=1 

J=1 



N 

Z (A jl A J3 ' 

N 

f A J1 A J3 } Z (A J2 A J3 + A J2 A J3 } 

<!>,/ 


J~1 

J=1 

J-l 



^Complex conjugate of the term (defined in Appendix B). 
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N 

^ ( b 3 A ji + b j A 3i > 

N 

£ ( b 3 A J2 + b J A 32> 

J=1 

N 

Z ( b J A J3 + Vfc> (20) 

J«1 

— 

Substituting the expressions for Aj 1# Aj 2 , Ajj, and bj 
into (20) results in three simultaneous equations which are solved for 
C Itj/ C D* and c Lg • 1116 algebraic expressions necessary to calculate 

the aerodynamic derivatives are presented in terms of measured response 
data. 



J-l 
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X <A 3l A J2 * A J1 A J2> * X [^- + “Vo Bln(|lS 2 - ^l 1 ] 


^ A J1 A J3 + A J1 A J3 * 


2q2s?B o . £ [ e o COB 0 2 " If Sin 0 l] 

A *•» ^ ^ - 


X ( Vj 3 * Vj3 )= '^T^ X K sln «*l] 


X (b 3 A Ji * b j A ji > ■ S! > S ) - %IVo '“Wi - 0i>] 


X (b ?J2 + V3 2 > ■ ° 

J*1 


X <b 3 A J3 + "jV ■ 2,S6 o X 


Kzz) z o 008 


Equation for pitching motion: 

Equation (19b) can be rewritten in the form, 


8 Jl C «e * «J2\ * 8 J3\ + *}k% * »j 
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where 

g j2 = * - HT 

“■“j 

«» - -.«f ••*** * ^ 

6jU = -qSc8 0 

hj = RlycP - Ke e )0 o e i?(2 ] J = 1,2,5, . . . N 


Applying a least-squares solution to equation (I9d) results 
in the following: 
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N 


1 

J=1 


N 


h (h J S J2 * Vfe* 

J=1 


N 


), (h J g J3 + h J g 33 ) 

J=l 


N 


), + h j«3^ 

fi 


(21) 


Substituting the expressions for g^, gj 2 > 6 jj» gj4» and hj into 
(21) results in four simultaneous equations which are solved for 0^, 
C^, C,^, and C^. The algebraic expressions necessary to calculate 

the aerodynamic derivatives are presented in terms of measured response 
data. 


» / ^ 2 \ 2 » 

h N" = (Vj i 

j=i j=i 


N 


L M 2 

J-i 


■( 


'qSff y 


2U / 


J=1 






20 o Zo (o . . , . 

+ — - — siu,0 2 - f^) 




O^O)- 
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N 


N 


>, M 

J=1 


2 _ 


(qSc)2 ) 

' 1 

J=1 


“I ♦(^J ■*■<* • V 


20 o z o ^ 




N 


'i |v| 2 - (^v M 

j-i 


N 


) (e jl g .i2 * - 2 (^) ) [ e o + ^ eln (« l 2.-« l l>| (“j) ! 

• i -J (mas* 


5 Jl B j2 6 jl 6 J2' " \ 2U 

J=1 J=1 


b (g jl 8 JJ + ‘jl^’ ' 1 ["Vo - *211^ 


J-l 


J=1 


N 


}_, + ‘ q 2U ^ } &®o Sln ^2 lL fo), 

j=l J=1 J 


M 


N_ 

), (g j2 g J3 + = 0 

J=1 



+ g J2^ : ) = - 


2q 2 S%5B c 


2U 


N _ 

£[■ 

j=i 


, o£z 0 , 
cu0 o sin 02 + cos 0^ 


■] 


O^CUj 


^ (g J3 g j4 

J=1 


n r 

+ ®j3®j 4 ) “ 2q 2 S 2 c 2 6 Q ^ e o cos 0 2 

J=i L 
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)_ (hJSj! + OjgJj) * 0 

J-l 


N 


N 


), (h J g J2 + h j«3 2 ) 

J=1 


) Bv* “ Kee^^o cos (^2 - <*i0 


OF 0)4 


N 


J=1 


(h j g J3 + h J g J3 ) 


N 

-2qgc ) 

j=i 


- k b. ) [« 


2 A <“0*0 . 
e o + — u“ sin (?>2 


-Cl)] 


a>=a)< 


N 


^ ( h j g j4 + h j6j4' “ -®9S»o ^ [ (I Y^ ‘ W 0 © cos <t] 


J=1 


J=1 


(JCPGD. 


The aerodynamic derivatives cure determined by solving 
equations (20) and (21). The data required are obtained from wind- 
tunnel tests. 

r 

At a test point (tunnel conditions fixed) the horizontal tail 
surface is oscillated at a known amplitude through a range of frequencies. 
At each discrete frequency co = cuj (j * 1,2,3, . . . N; N > 4) the model 
notion is monitored to determine z Q , e o , and $ 2 * Cable tensions 

are monitored at each test point to calculate mount stiffnesses from the 
equations presented in Appendix A. The data collected at N different 
frequencies are used to evaluate the algebraic terms which comprise 
equations (20 ) and (21). Bach equation is then solved for the associated 
aerodynamic derivatives. 
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5.2 Error Analysis 
5.2.1 Error Analysis Computation 

An error analysis is provided to determine the accuracy with which 
the model response z D , 9 0 , fa, and fa needs to be measured in order 
to obtain meaningful results. A numerical error analysis is presented to 
determine the error in the calculated aerodynamic derivatives due to an 
error in measuring model response. Since no experimental data are avail- 
able, a set of amplitudes and phase angles is obtained by solving the 
equations of motion (eqs. (17a) and (17b)) for an assumed tunnel test 
point and representative values of th aerodynamic derivatives This 
analysis assumes that the governing equations of motion are correct as 
presented in equations (17a) and (17b). All parameters other than model 
response are assumed to be correct. 

Errors are now introduced into the response data, and the aero- 
dynamic derivatives based on the new response are calculated using 
equations (20) and (21). The new aerodynamic derivatives c~n be 
expressed in a Taylor series in the following manner: 


& 


- r J. 8c . . 8c 

■ -0 dez e z a €e 


Ae 


0 8e, 


01 


Ae 


0i 


8c 

8e 


Ae 


fa 


02 


( 22 ) 


1 The aerodynamic derivatives selected are those of a large 
multi Jet cargo airplane. 


where 
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= numerical value of aerodynamic derivative 
with error 

= numerical value of aerodynamic derivative 
without error 

= incremental error in measurement 

= rate of change of the aerodynamic derivative 
with respect to the erro > * 

Appendix C describes the computer programs used in this analysis. 

The model properties, test conditions, and aerodynamic derivatives 
assumed for the numerical error analysis are as follows: 

U = 500 ft/sec 

q =100 lb/ft 2 

S = 10 ft 2 
m =2.0 slugs 
c = 1.0 ft 
Iy = 5 .0 Slug -ft ^ 

6 q = 0.00698 rad 
= 20.0 lb/ft 
Kqq = 500 ft -lb/rad 


a) = 1,2,3, 


• • • 


30 rad /sec 


Cj^ = 5.00/rad 
C D = 0.02 

= -1.00/rad 
Cj^ = -4 .oo/rad 
= -15.0/rad 
Cj^ = -0.4o/rad 
= 1.2/rad 


vfo 


C €=0 


Ac 


dC 

8c 


4o 


These parameters are substituted into equations (C-l) and (C-2) to 
evaluate the model response Z Q , 0 Q , and as a function of the 

tail frequency gd. The model response for this example is presented in 
Figure 4 . The peak responses occur at the damped natural frequencies of 
the system. 

Errors are now introduced into each term of the response data, and 
the aerodynamic derivatives with response errors are calculated using 
equations (20) and (21). Appendix C presents equations (20 ) and (21) as 
equations (C-3) and (C-4), where 


N 


B(l,l) = ^ | Ajl | : 
>1 


B ( 1 >2) = g ^ ( A jl A j2 + A J1 A J2 ^ 

J=1 


N 


B(l,3)=i £ (A^^vy, etc. 

J=1 


N 


AU, 1 ) - ^ |gji| 


N 


A(l ' 2) - i 1 (8 3l g J2 + V&* 


J=1 


A(l ,3 ) = | ^ etc. 



Figure 4.- Longitudinal response as a function of tail frequency 
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Based on equation (22), a numerical comparison of the aerodynamic 

derivatives, with and without response errors, is used to calculate — 

de 

and the per cent error. In general, equation (22) is written 


Vo = *5^2 *1^5 ♦ 




Se, 


Ae 


m 


m 


Considering each of the error functions separately (Ae^ ■/ 0 } 

Ae 2 = Aej . . . = Ae m = 0) results in the following expressions for the 
rate of change of the aerodynamic derivatives with respect to the 
response errors: 


dC _ C ei/0 ~ C €=0 dc _ Cg 2^0 " C e =0 

de^ Ae^ dc 2 Ae 2 

We can also define 


( 2 ?) 


Per cent error = Z.? € . z9 x ioo (24) 

C e=0 

A numerical calculation of the aerodynamic derivatives with and 
without response errors is presented in Table I. Errors ranging in 
amplitude from 0 to 5 per cent and in-phase angle of 1°, 2°, and 3° are 
presented. The derivatives are obtained numerically by considering an 
error in one response function while the others are held constant. 

Tables II and III present the calculated data in terms of the 
rate of change of the aerodynamic derivatives with respect to response 
error and the per cent error in the aerodynamic derivatives, respectively. 


TABLE I.- CALCULATED AERODYNAMIC DERIVATIVES WITH RESPONSE ERROR 


rad 

C D 

5.00 

0.02 

5.05 

-0.03 

5.10 

-O.08 

5.15 

-0.13 

5.20 

-0.l8 

5.25 

-0.23 


w € e» 

fs error 


4.95 

0.07 

4.90 

0.12 

4.85 

0.17 

4.81 

0.21 

4.76 

0.26 


Cn a 

c md. 



+ c me 

rad 

rad 

rad 

rad 

rad 

-1.00 

-4.00 

-15 .00 

1.20 

-19.00 

-1.00 

-3*75 

-15.05 

1.197 

-18.80 

-1.00 

-3.42 

-15.17 

1.194 

-18.59 

-0.99 

-3.01 

-15.37 

1.191 

-18.38 

-0.99 

-2.54 

-15.63 

1.185 

-18.17 

-0.99 

-2.00 

-15.97 

1.182 

-17.97 


-1.00 

-4.17 

-15.03 

-1.00 

-4.26 

-15.14 

-1.00 

-4.27 

-15.32 

-1.00 

-4.20 

-15.58 

-1.00 

—4 .06 

-15.91 


■9.20 

■9.4o 

•9.59 

.9.78 

■9.97 


5.00 

0.09 

5.00 

0.17 

5.00 

0.24 


-0.98 

0.96 

-19.89 

-O.96 

5.78 

-24.61 

-0.94 

10.46 

l 

-29.14 


5.00 

-O.06 

5.00 

-0.l4 

5.00 

-0.22 


-1.04 

-7.53 

-1.09 

— 18 .86 

-1.13 

-26.25 
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TABLE II.- RATE OF CHANGE OF AERODYNAMIC DERIVATIVES WITH RESPONSE ERROR 
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TABLE III.- PER CERT ERROR OF AERODYNAMIC DERIVATIVES 


% error 

C Iu 

C D 

% 


C ®a 

C mg 


^ + 

1 

1.0 

-250 

1.0 

0 

-6.3 

1 

-0.25 

-1.1 

2 

2.0 

-500 

2.0 

0 

-13.5 

■ 

-0.50 

-2.2 

3 

3.0 

-750 

3.0 

•1.0 

-24.8 

2.5 

-0.75 

-3.3 

4 

4.0 

-1000 

4.0 

-1.0 

-36.5 

4.2 

-1.25 

-4.4 

5 

5.0 

-1250 

5.0 

-1.0 

-50.0 

6.7 

-1.50 

-5.4 

7 > error 









i 

-1.0 

250 

0 

0 

4.3 

0.2 

1.3 

!.l | 

2 

-2.0 

500 

0 

0 

6.5 

0.9 

2.5 

2.2 

3 

-3.0 

75' 

0 

0 

6.8 

2.1 

3.8 

3.1 

4 

-3.6 

900 

0 

0 

5.0 

3.9 

5.0 

4.1 

5 

-4.8 

1200 

0 

0 

2.0 

6.1 

6.3 

5.1 

deg 









1 

0 

350 



-124.0 

33.0 

-0.50 

-0.37 

2 

0 

750 

-36.0 

-4.0 

-245.0 

64.o 

-0.75 

-0.90 

3 

0 

1100 

-54.3 

-6.0 

- 362 ,0 

94.0 

- 1.00 

-1.70 

Ae^p, 

deg 

■1 

u 







■■1 

■ 



4.0 

88.0 

-24.0 

0.75 

*0 .2x 

2 

0 

Q| 


9.0 

374.0 

-100.0 

1.50 

-0.95 

3 

EH 

91 

i/ 

13.0 

556.0 

-151.0 | 

2.25 

-2 ,10 











These Tables sure obtained by applying equations (?3) and (2k) to the 
results presented in Table I. Error intervals are included to establish 
the linearity of the functions. Table II is present’d so that equa- 
tion (22) can be used to establish the effect of errors in several of 
the variables on each aerodynamic derivative. 

3.2.2 Results and Discussion 

The results given are for a particular example; therefore, they 
should not be generalized. Since the dynamics of the problem are 
dependent on the system parameters and the test point, it is necessary 
that each case be analyzed separately. 

Referring to Tables I and III, notice that Cj^ appears to behave 
favorably in the dynamic analysis. That is, Cj^ is not significantly 
affected by errors. This term is one of the more important derivatives 
and is normally measured statically. Numerical values of this derivative, 
with errors in response amplitudes, tend to give linear results with error 
in the derivative on the same order of magnitude as the error in the 
response. Table I shows that errors in the phase angles ^ and $ 0 
do not affect this derivative. 

(2) C D 

Referring to Table I, large errors in the drag coefficient are 
evident for all values of response error. The reason for this is 
apparent from equation (15a). The drag coefficient appears in combination 
with in the form (Cj^ + C D ) . For most configurations Cj^ > > Cp, 

and as our results show, the drag coefficient cannot be separated when a 
response error is introduced. 
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Although the drag coefficient, then, should not be measured using 
the dynamic technique, it is still possible to obtain Cp statically on 
the mount system. If the static equations of equilibrium are written in 
the x direction (refer to Fig. 3), the front and rear cables differ in 
tension as a function of geometry and model drag in the following 
approximate manner: 


hence 


2Tp cos = 2T R cos 0 R + qSC^ 


C D * 


2(Tp cos 0p - Tr cos 0 r ) 
qS 


Load cells on the front and rear cables provide Tp and Tp. 

(3) ^ 

Referring to Table III, note that amplitude errors affect the 
magnitude of Cpg slightly, ranging from zero to about 5 per cent. 
Errors in phase effect this term appreciably. Errors of 3° in either 
^ or $2 result in magnitude errors of over 50 per cent . 

The reason for this result is that for conventional airplane 
configurations the tail lift is almost negligible. Hie primary forcing 
function obtained from the tail is the aerodynamic moment M^. This 
moment is physically generated by the tail lift acting a distance 

I>t rearward of the model center of gravity. (i^. is the distance 
between the tail aerodynamic center and the aircraft center of gravity . ) 
Etkin (Ref. l) shows that the aerodynamic moment can be approximated by 
the equation - -L^Z^ . 
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Hence 


qScCngS^ -LtqSCj^B 


or 


QngC 


For our example 1^. = 3.0 Ft since "c = 1.0 ft, = 1.2/rad, and 

-0.4/rad. We can therefore assume that Cj^ can be determined 
from the values obtained for • 

( 4 ) 

Referring to Table III, notice that behaves favorably in 

the dynamic analysis. Amplitude errors affect the magnitude of 
slightly, ranging from near zero to about 6 per cent. Error in phase 
angle affects even less. 

(5) 

Another of the important derivatives, which is normally measured 
statically, appears from Tables I and III to behave quite well in the 
dynamic analysis. The effect of amplitude errors on is quite 

small. Errors in phase angle affect this term the most, resulting in a 
13 per cent error due to an error of 3° in 

(6) CmQ > (^ny + ^mg) 

Tables I and III show that errors in both amplitude and phase 
measurements result in extremely large errors when C,^ and Cj^ are 
calculated separately. For the test condition analyzed, errors of over 
500 per cent are evident for small errors in phase. 
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Notice, however, that the sum of these terms (C,^ + C^) behaves 
quite well for both amplitude and phase errors. Numerical results 
presented in Table III show that the sum (Cm^ + C^) varies almost linearly 
with amplitude errors while phase affects the results a maximum of about 
2 per cent. In view of these results, it appears that the sum of 
(Cjd^ + Cjd^ ) can be measured accurately in the dynamic analysis, but 
unique solutions for each of the derivatives cannot be determined. The 
following section will farther interpret this problem. 

5.2.3 Comparison of Free-Flight and Two-Cable-Mount Equations of Motion 
In order to fully understand the problems which arise in deter- 
mining unique solutions for each of the aerodynamic derivatives, a 
comparison of free-f light and mount equations is presented. Greenburg 
(Ref. 2) shows that for the free-flight case a linear dependency exists 
in the dynamic response, so that the aerodynamic derivatives in the pitch 
equation cannot be solved uniquely. Neglecting the drag coefficient Cp 
the free-flight lift equation (eq. (15a) , K? 7 = 0) can be written as 

follows : 

♦ 

mz + qSC T £ + 

■4x U 

Since a = 9 + equation (25) can be written as 


qSC^Q — 


(25) 


mUa - mU9 + I^a + I^& = 0 


(26) 



where 
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h. - qisc^ 

Lg » qSCj^ 

Similarly the free-flight pitch equation becomes 


lye - 3§sf(p + cl.) 0 - qScc_ e - c .z - — c i 
™ 2U % “0 H ®a 2U 2 “a U % 


or 


Where 


I Y® - *t a ' " “e® = M 5 B 




% = 


qSc 2 

2 r 

2U 


a 


qSc 2 

M 0 = ~W C m • 


M& = qscCL 


Solving equation (26) for e results in 


. I^a + Lg5 

0 s + a 

mU 


= qScC^S 


(27) 
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Substituting this expression into equation (27) results in the following: 

V + + - 8 (% t! ^)- 0 < 23 > 

Due to the relation expressed in equation (26), equation (28) can be 
solved for the following combination of derivatives: 




mU 


MqLb 


Mr + 

5 mu 




Greenburg (Ref. 2) states that for most conventional conf igurat ions this 
indeterminacy affects only the separation of the damping derivatives 
(C^ + ) because they are of the same order of magnitude. For most 


configurations the terms + 


Meifa 

mU 


and Mg + 


mU 


are approximately 


^ and Mg, respectively, to within 5 per cent. Only if Mq or 
can be determined separately by some other method can the damping 
derivatives be separated uniquely. 

The equations equivalent to (26) and (27) on the mount system can 
be written as follows: 


mUd - mU0 + I^a + Lg8 + K^z = 0 


(29) 
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Iy'e - %a - J^a - + KqqO = Mge ( 30 ) 

Solving equation (29) for 0 results in 

. Iu. a + l 5 6 + K ZZ Z 

8 = + a 

mU 

Substituting this expression for 0 into equation (30) results in 


I Y 0 + KQQ0 - am + 


d ^ + M e> + ^ 


Note the similarity in the aerodynamic ^erms comprising equations (28) 

and (31) • If the dependency derived in the free-f light case were to 

exist in the mount system analysis, the case where exact data are used 

would also break down. It has already been demonstrated in Table I that 

the mount analysis yields each of the aerodynamic derivatives uniquely 

M/SK77Z 

when exact response is used. It is the added expression — „ which 

mu 

allows us to solve for uniquely and, therefore, solve for each of 

the aerodynamic derivatives. Hence the mount system restraint K^ z 
allows both the vertical translation and pitch equation to be solved 


separately for the aerodynamic derivatives C^, C^, C^, C^ & , C L& , 

and Cj^ . 


The analysis fails to separate the terms 


Si + Sj 


when an error 


is introduced into the response data due to the magnitude of the term 
M/K^z 

- „ . For example: 

mu 
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3 Sc 2 


=101 


a?z 0 ei ^ 1 -2 
t qSc results in 


Combine — C m .K !z W *ith ^ ^ 

Yjzjl o 

approximately -jj- compared to or. For this example = 20 lb/ft, 

m = 2.0 slugs, and a> = 1,2,3, ... 30 rad/sec. Hence for a large 
range of a>, the expression defining C^. separately is small compared 
to the expression defining (Cj^ + Cj^), and the introduction of an 
error into the response eliminates any accuracy in separating the 


two terms . 


CHAPTER IV 


EXPERIMENTAL TECHNIQUE FOR MEASURING AILERON EFFECTIVENESS 
AND DAMPING-IN-ROLL STABILITY DERIVATIVES 

4 . 1 Int roduct ion 

The lateral equations of motion of the two -cable -mount system are 
given by equations (l4c), (l4d), and (l4e). The problem of determining 
each of the aerodynamic derivatives has greatly increased over that of 
the longitudinal case due to the added degree of freedom and its asso- 
ciated derivatives. 

The general approach to measuring the three-degree-of -freedom 
lateral derivatives would be the same as discussed in the previous 
chapter. However, Etkin (Ref. 1) states that for many conventional 
airplane configurations, the roll equation in free flight can be simpli- 
fied and treated as a single degree of freedom. If we assume it is 
possible to force the model to behave as a single degree of freedom in 
the wind tunnel, equation (l4d) could be simplified as follows: 

- (ilr % - 

Letting L^ = qSbCj 6^ results in 

" (^" C l p " + = qSbC Z 6 8 A (52) 
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The aerodynamic derivatives Cj and Cj appearing in 

6 p 

equation (}?) are referred to as the aileron effectiveness and damping- 

in-roll stability derivatives, respectively. A brief insight into the 

physical nature of the derivatives C 7 and C, is necessary. 

l b 4 p 

The primary function of the ailerons is to produce a rolling 
motion of the aircraft. As the tern implies, aileron effectiveness is 
a measure of the performance of the ailerons in producing this motion. 
Physically, a differential deflection of the ailerons creates an incre- 
mental change in the lift on each wing. Since this lift ie in opposite 
directions on each wing, a rolling moment is produced. However, when 
dealing with an elastic airplane, a deflection of the ailerons also 
produces a twist of the wings. This Induced twist changes the wing lift 
in the opposite direction to the lift due to the aileron deflection. 
Hence, the moment generated by a control input is the difference between 
the moment produced by the ailerons and that induced by wing twist. An 
elastic airplane may even experience a phenomenon referred to in the 
literature as "aileron reversal," if the moment generated by the wing 
twist is larger than that produced by the imposed aileron deflection. 

is referred to as the damplng-ln-roll derivative. In most 
configurations only the wing contributes significantly to :Ms de*' va- 
tive. A rolling moment is generated opposing the rolling motion of an 
airplane due to the spanwlse angle of attack produced by a roll rate. 

The angle of attack varies linearly across the wing, from a value of 
|jj at the right wing tip to - ~ at the left wing tip (+# right 
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wing down) . The angle-of-attack distribution along both wings generates 
a moment which resists the rolling motion of the airplane. 

It will be shown in this chapter that through the proper selection 
of cable-mount parameters, the roll response of the model to a sinusoidal 
oscillation of the ailerons can be approximated by a single degree of 
freedom. A parametric study of the model and its mount system is pre- 
sented to establish a mount configuration which permits the roll response 
to be approximated by a single degree of freedom. 

Once the single-degree-of -freedom system approach tc the problem 
is established, the equation of motion is solved for the aerodynamic 
derivatives and C^, based on the dynamic response technique 

presented in Chapter III. In this case the ailerons are sinusoidally 
oscillated through a known frequency range, and the steady -state model 
response Is monitored. 

In order to verify the dynamic technique, experimental result 5 
are presented for an aeroelastically scaled model of a high-speed Jet 

transport that was tested in the wind tunnel to determine Cj and C j . 

8 P 

IXie to the static nature of C it is possible to measure this deriva- 
tive without resorting to a dynamic approach so that a static measure of 



can be found for comparison with the dynamic results . 


The mount 


system used during the static wind-tunnel tests is also unique and will 


be described. Since Cr is a dynamic derivative, no experimental 

P 

results are available for comparison purposes. A simplified error 


analysis to determine the accuracy of the dynamic technique is also 


presented. 
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4.2 Analytical Aspects 
4 . 2.1 Equations of Motion 

The lateral equations of motion on the two-cable-mount system are 
given by equations (l4-c), (l4d), and (l4-e). In order to measure the 
aileron effectiveness and damping -in-roll stability derivatives, it is 
necessary to show that the model response to a sinusoidal oscillation of 
the ailerons can be simplified to a single degree of freedom. If we 
assume that an oscillation of the ailerons only generates a forcing 
function in roll, then the right-hand side of the lateral equations of 
motion can be -itten, as a function of time, in the following manner: 

L(t) = qSbC 2s 6 A e ia * 

N(t) =0 (35) 

Y(t) = 0 

Assuming sinusoidal motion, the steady-state response becomes : 

«(t) - 0 o « 1( “* <a l) *<t) - y(t) - 

j{t) = t(t) - la» 0 e 1(alt ^' y(t) = 

jf(t) - *(t) . y(t) - -^y^^ 

The phase angles relating model response to aileron deflection are 

o»l , 02 , and • 

Substituting equation (33) and the steady-state response into the 
lateral equations of notion results in the following form of the 


equations: 
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Roll: 

^o e 1 j" I x a ^ + ^ " i(a i^T Cj p “ 2 ^ I x a #0)jj + +o e + ^ Sbc Zp 

+ %- iUiS W C J + y o eia3 [vY - i<u lT %] = qSbC % 8 A 

(54) 


Yaw: 


0o eiai haft + m - ia > ^ c npJ + t D e la2 -I Z a? + qSbC^ 




ia3 ( 3 fu“ ° n r " 2 ^ I Z u wjJ + yoe^^Y - i«> C np J = 0 


(35) 


Si ~e translation: 


^ -^-^W c y p ] + *o ela2 [ qSbc y B + K » * 110 W c y r ] 


y Q e 


lai 


[? ( % ' 


-TOO? + Kyy - ico 77k Cy_ - C D ) - 2^ymayy 


Y 


= 0 


( 36 ) 
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The simplification of equation (54) as a single degree of freedom 

yields, 


s iai 



-Ixu? + - iu>(— - C 7 - 1 - qSbC 3 & 5 a ( 57 ) 


)}■ 


It has been stated previously that in flight the roll equation can be 
analyzed in a form similar to equation (57). Since the cable restraints 
are a function of mount system geometry and other test parameters, it is 
necessary to compare calculated results based on both single degree and 
three-degree-of- freedom solutions. Assuming approximate values of the 
aerodynamic derivatives, equations (54 ) , ( 55 ) , and ( 56 ) are solved 
simultaneously for the roll response as a function of the forcing fre- 
quency oj. These results are then compared with those for the "single 
degree" equation. A parametric study is presented for a numerical 
example with physical and mount properties quite similar to the model 
of the high-speed jet transport tested in the wind tunnel. 

4.2.2 Parametric Study of the Model and Its Mount System 

Let Mach number = O. 89 , q = 225 psf, hence U = 470 ft/sec. 

The model physical properties, mount configuration, and assumed aero- 
dynamic derivatives^" are as follows: 

^The aerodynamic derivatives given are estimated through informa- 
tion supplied by the aircraft manufacturer. 
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2.18 slugs 

% 

= -0 .4oi/rad 

5.25 slug-ft 2 

% 

= 0.078/rad 

7-50 slug-ft 2 

% 

= - 0 . 062 /rad 

0 

% 

= 0.422/rad 

8.94 ft 2 

C "r 

= -0.124 /rad 

8.46 ft 

% 

= 0.117/rad 

0.55 ft 

% 

= 0.105/rad 

0.50 ft 

S r 

= 0.005l/rad 

3 R = 20.0 deg 

s 

= - 0 . 725 /rad 

Lp = 20.0 ft 

C D 

= 0.02 

l6l lb 


= 0.035 

l4o lb 



Cy = Cy = 0.05 




The parametric study will ^eal with the design of a mount con- 
figuration based on the fore and aft pulley separation distances a 
and e (Fig. 5), which permits the roll response to be approximated by 
a single degree of freedom. The cable restraints for this example are 
calculated from the mount restraint equations presented in Chapter II. 
Calculated values of the mount restraints Kyy and are 28.5 lb/ft 

and 70.6 ft-lb/rad, respectively, since these restraints are not functions 
of either a or e. Values of the restraints Kyy and Kyy, as a 
function of the parameters a and e, are presented in Figure 5. Due 
to the mount symmetry Ky^ = K^y = Ky^ = K^y = 0. 
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Using the physical and aerodynamic properties of the example 

given, equations (3M, (35), and (36) are solved simultaneously for 

the roll response as a function of the aileron frequency co. Figure 6 

presents results obtained for three sets of a and e values and 

compares three -degree -of -freedom to single-degree-of -freedom solutions. 

Since no value of C? c is assumed, the results are presented in terms 

0 

of the magnitude — — — and versus aileron forcing frequency to. 

C Zg °A 

These results should be valid for any nonzero C, . 

-6 

All three sets of data presented in Figure 6 agree well with the 
single-degree-of -freedom calculations for frequencies above approximately 
8 rad/sec. Below this frequency, a peak response around 4-l/2 rad/sec is 
apparent for a = 0, e = 1.75 ft, and a = 0, e = 1.0 ft. For the 
case a = 1.2 ft, e = 1.0 ft, this peak response is much less pronounced 
and the calculated data compare favorably with the single degree analysis 
throughout the frequency range given. 

for the three sets of a and e values (Fig. 5) are as 

follows : 

- ft -lb/ft 

a = 0, e = 1.75 26.6 

a = 0, e = 1.00 l4.8 

a = 1.2, e = 1.0 0 

Comparing the trend in with the results presented in Figure 6, we 

notice that as the magnitude of diminishes, the single degree and 

three-degree-of- freedom analysis agree more favorably. 
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These results are more apparent if we look at the homogeneous 
solution of the equations of motion to determine the dynamic behavior 
of the system. Let L(t) = N(t) = Y(t) = 0 . Assume 0(t) = 0 o e^, 
t(t) = t 0 e At , and y(t) = y 0 e A * ; where A = rj + ire. Hence 

#(t) = A0 Q e At t(t) = At c e At y(t) = Ay c e At 

#(t) = A 2 $ D e At t(t) = A 2 >|r 0 e At y(t) = A 2 y Q e At 

Substituting these terms into equations (3*0, (35), and (36) results in 
three homogeneous algebraic equations in the unknowns 0 O , y c , 

and containing A . 

Therefore 

tfo&ll* 2 + A 12 A + A 13D + ^o[ b ll A2 + b 12 A + ^13] + y o[ C 12 A + C 13 D = 0 
^oC A 21 A2 + A 22 A + k 2^] + ♦o[ b fel A2 + b 22 A + b 23D + y o &22 A + C 2y\ = 0 
^o[ A 32 A + a 333 + + o[ b 32 A + b 33j + y o &31 a2 + C 32 A + C 33D = 0 


( 38 ) 



qSb 
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Getting the determinant of the coefficients of equation (38) equal to 
zero provides the condition for determining the values of A. This 
determinant is known as the "stability determinant," and its expansion 
results in the characteristic equation of the dynamic system. For this 
example the characteristic equation is of the form 

AA 6 + BA 5 + CA^ + DA 5 + EA 2 + FA + 0 = 0 (39) 

The six roots of this equation establish the dynamic characteristics of 
the system. For the example given the solutions of equation (39) are as 
follows : 


a = 0, e = 1.75 

A x - 

-1.50 


II 

£ 

-11.1 


A 5 ,Au = -0.16 ± i 4.1+7 
A^,A6 = -2.28 ± i 18.79 

a = 0, e = 1.0 A x = -1.407 

A 2 = -11.16 

A,,A 4 = -0.152 ± i 4.15 
?yA 6 = -2.22 ± i 17.9 

a • 1.2, e * 1.0 = -1.17 

Ag = -11.17 

= -0.56 ± i 3-6 
A^,Ag - -2.27 ± i 19.2 
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Figure 7 presents these results plotted in the complex plane. Since the 
complex roots appear in conjugate pairs, only the upper half of the 
complex plane is presented. The radial distance from the origin to a 
complex root is the undamped natural frequency o^. The angle between 
the radial vector and the ioo axis is equal to sin~l £ n , where £ n 
is the damping ratio of the model in flight relative to critical damping. 

Note in Figure 7 that a lightly damped mode exists at 4.47, 4.15, 
and 3-6 rad/sec for the three cases under consideration. Reed and Abbott 
(Ref. 4) refer to this as the side translation mount mode. (This can be 
verified by substituting the roots A^, A^ into equation ( 38 ) and 
solving for the characteristic mode shape.) If we assume the forced 
model response at the side translation frequency decreases with increased 
damping in this mode, the case where a = 1.2, e = 1.0 would exhibit 
the least response. Figure 6 confirms this fact since the response at 
the side translation frequency for a = 1.2, e = 1.0 is much less 
pronounced. The cable configuration which allows single degree approxi- 
mation seems to be one in hich the mount restraint K^y = 0 since this 
increases the damping in the side translation mode. Based on this 
analysis the model pulley locations were fixed at a = 1.2 ft, 
e = 1.0 ft, and the model response to a sinusoidal oscillation of the 
ailerons was assumed to behave as a single degree of freedom in roll . 
4.2.3 Single-Degree-of -Freedom Solution 

Once it has been assumed that the model behaves essentially as a 
single degree of freedom in roll, equation (37) can be written in the 


form 
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Figure 7 .- Dynamic characteristics of lateral modes 
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^(^ir % - p Vx^f)^o eiai + ^ sbc z 6 6 A = ^- I x a)2 + K ^o eiai 

( 40 ) 

One further simplification is desirable in order to measure the aileron 
derivative Ci and the damping-in-roll derivative C? . Looking at the 

8 V p 

magnitude of the terms — — and i n equation (4o), we 

pl2 P 

notice that (Cj « 0.3 from estimated data). 

Therefore, equation (4o) can be written approximately as 


10 ) 


qSb £ 

2U 


^2 aic i 


«S OCi^A - (-Ij^ 


iai 


(41) 


Equation (4l) can be solved for and using a least- 

p 8 

squares solution, by measuring the dynamic amplitude <$ Q and the phase 

angle at N discrete frequencies. Assuming and are 

8 p 

independent of frequency, equation (4l) becomes 


A jl C i T 


A J2 C l* 


= b . 


■where 




\J1 " 2U 


OJ=6D4 


Aj 2 = qSb8 A 


A J5 » [(.1^ <■ V0 o e la i] j ■ 


N 
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From Appendix B 


N 

N 

2 Z M 2 

C\J 

<p 

t— 1 

J=1 

j=i 

N 

N 

^ A Jl A j2 + A jl A j2^ 

OJ 

cvj 

< 

OJ 

J-l 

j-i 


+ A ,11 A j2 ) 




1 (A ji b j* A ji b j> 

J-l 

)_, ^ A j2 b j +A j2 b j) 

J=1 


Hence 


r oX 2 N_ 

ss£) y w f 

2U J /_, v ^o''oj=cuj 

J-l 


q 2 S 2 b36 A 


N 


2U 


) (cu0 sin a-, ) 

/_, ° l'cD=GOj 

J=1 


q 2 S 2 b?5 A f , J 0 

2U K> Sina l^a„ ^ Sb M N 


J-l 


0 


N 

-qSb6 A ^ [(ijftiP - K00)(0 q cos 

J-l J 

(42) 


s 
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4 .3 Experimental Results 

A l/l9-size aeroelastically scaled model of a large, subsonic, 
multijet cargo airplane was tested in the NASA-Langley transonic 
dynamics tunnel. The wind tunnel uses Freon-12 gas as the testing 
medium. The speed of sound in Freon-12 is approximately half that of 
air. Since many aeroelastic phenomenon are functions of Mach number, 
this allows testing to higher Mach numbers at reduced dynamic pressures . 

Since can be measured either statically or dynamically, 

test results employing both of these methods will be presented. These 
results will then be compared to establish an evaluation of the dynamic 
approach . The mount system used during the static tests is also unique 
and will be briefly described. 

4.3.1 Static Wind-Tunnel Tests 

Once the model is designed and constructed, static testing is 
fairly straightforward. Grosser (Ref. 3) presents a testing technique 
where the model is supported in the wind tunnel on a sting-pylon-spring 
mount system. A photograph of the model on this mount support is 
presented in Figure 8. 

The support consists of a rigid sting which is attached to the 
tunnel splitter plate, a pair of pylons which extend from the sting to 
within the model fuselage contour, and a set of springs which connect 
the fuselage spar to the pylons. The mount allows the model six degrees 
of limited freedom . The model is restrained from excessive motion by 
stops located on the front and rear pylons . Rolling moments generated 
by a deflection of the ailerons are transmitted by a pair of push-pull 



Photograph eh model ot st ing-pyio'. -jpring mourl 
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rods extending from the fuselage spar to a strain-gaged beam located in 
the sting. Roll control for trim, independent of model control surfaces, 
is provided by a remotely controlled rotation of the entire roll 
mechanism within the sting. 

At selected test points during a run, the ailerons are remotely 
deflected. Rolling moments and aileron displacements, measured from the 
strain-gaged beam and position indicators on the ailerons, respectively, 
are monitored on direct writing recorders. Since rolling moment is 


defined as qSbC 7 5a# C 7 can be determined from the measured data. 
L b H L b 

Static values of C 7 , determined at the same test conditions as the 


dynamic tests, will be presented under the comparison of experimental 


results . 


The static wind-tunnel tests were originally designed to establish 
the aileron reversal boundary of the model. For purely informative 
purposes the results of this test are presented in Figure 9- Curves 
representing the measured reversal boundary (C 7 = 0) and 25 in-lb of 

rolling moment/degree of aileron deflection are given as a function of 
model dynamic pressure versus Mach number. The two reversal points 
measured statically on the cable-mount system will be discussed later. 

The dashed line representing the estimated boundary is based on positive 
aileron effectiveness measurements approaching reversal. The reversal 
boundary could not be attained due to excessive buffeting loads experi- 


enced in this region. 



Figure 9 *- Aileron effectiveness boundaries 
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4 dynamic Wind-Tunnel Tests 

In order to verify the dynamic technique, the model used in the 

previously mentioned static tests was tested dynamically to determine 

the derivatives C, and C, . A photograph of the model on the two- 
£ & £ p 

cable mount is presented in Figure 10. 

The aileron drive mechanism on the model had to be modified in 
order to provide a sinusoidal fixed amplitude deflection. The aileron 
drive was modified to consist of a pair of push-pull rods extending 
from the aileron pivot to the wing attachment structure, a rotating cam, 
and a variable speed ac motor. The aileron frequency is altered remotely 
by varying the voltage to the drive motor. A finite torque was required 
to overcome friction in the drive system, and therefore limited it to a 
minimum sustained frequency of about 0 .5 cps . The maximum frequency was 
around 4 .0 cps . The aileron amplitude could not be altered during a run 
since this was set by the eccentric attachment of the push-pull rods on 
the cam. The drive mechanism was designed so that the ailerons would 
return to zero deflection after the oscillation. The ailerons could also 
be statically deflected up to the maximum preset dynamic amplitude. Roll 
control was provided remotely by a pair of spring loaded spoiler panels 
located inboard on each wing. These panels were opened and closed by 
means of a torque tube arrangement driven by an electric motor. Longi- 
tudinal control was provided by a remotely controlled horizontal 
stabilizer. 

Onboard instrumentation included a miniature rate gyro and a 
strain-gage position indicator to measure roll rate 0 and aileron 



ograph ot* mod*.- 1 or. two- w it 1 * m ant 
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displacement 5^, respectively. A servo accelerometer was included to 
measure static roll angles and provide a level roll reference. Signals 
from these instruments were displayed on a direct writing recorder for 
visual monitoring and simultaneously recorded on magnetic tape for 
further analysis. 

The dynamic amplitude is determined by integrating the roll rate 
signal. For sinusoidal motion this integration results in |0 O | = 
Ideally, a direct comparison between the traces $ and as a 

function of time is all that is required to determine the phase angle 
• Since the model is essentially free and subjected to tunnel turbu- 
lence, signal noise, etc., this procedure is quite time consuming and 
subject to added errors. In lieu of this, an electronic sin-cosine 
resolver was utilized to determine a^. This instrument is designed to 
electronically evaluate the phase angle between two known signals. 

Once the dynamic response of the system is known, a value of the 

mount restraint is required before solving equation (42) for the 

derivatives C, and C, . This restraint in calculated based on the 
l b l p 

equation presented in Chapter II. Since is a function of front 

and rear cable tensions (geometric properties are known), miniature load 
cells are installed in these cables and at each test point these tensions 
are recorded. 

The geometric and physical properties of the model and mount are 
as follows: 1 



1 A response analysis with measured model properties was run to 
establish that the model behaves essentially as a single degree of freedom. 
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I x -- 2.16 slug -ft 2 
Weight = 47.9 lb 
S =8.94 ft 2 
b = 8.46 ft 

= 0.105 rad 
Lp = Lp = 25 .0 ft 
3 F = 3 r = 20 deg 
a = 1.2 ft 
e = 1.0 ft 
h =0.57 ft 
d = 0.59 ft 

Tests were run at Mach numbers 0.675 and 0.75 over a range of 
dynamic pressures from 115 psf to the reversal boundary. At each 
dynamic pressure, tunnel conditions were held constant while the ailerons 
were oscillated over a range of frequencies. At each discrete frequency, 
roll rate, aileron deflection, cable tensions, and the tunnel parameters 
were recorded. 

Theoretically, the model should have no roll response at the 
reversal boundary. Such is not the case, since this assumption is 
based on a single-degree-of -freedom analysis. For the areas where the 
ailerons are quite effective is a function of Mach number and 

dynamic pressure), the assumption that N(t) = Y(t) =0 is realiztic. 

As the ailerons become less effective 60 does the forcing function L(t). 
In the region near reversal the servo accelerometer was used. 



ai 


The reversal dynamic pressure was estimated at the point where no roll 
was evident for a static deflection of the ailerons. The dynamic pressure 
was then increased slightly above this point, and reversal was noted both 
on the accelerometer and visually. 

Table IV presents the measured model response data as a function 
of Mach number and dynamic pressure over a range of aileron forcing 
functions. The dynamic amplitude 0 O , the dynamic phase angle a^, 
(sin-cosine resolver), the cable-mount tensions, and the tunnel test 
parameters at each frequency are substituted into equation ( 42 ) to deter- 
mine the aerodynamic derivatives. Figure 11 presents the calculated 
aerodynamic derivatives as a function of Mach number and dynamic pressure. 

4.3 .2.1 Error Analysis .- In order to establish the effect of 
measured response errors on the accuracy of the calculated derivatives, 
a simpl numerical error analysis is included. Test data obtained at 
Mach number 0.675 are used in the analysis. At each of the test dynamic 
pressures a response error was introduced in the following form: 


*Wo = ^=0 + e 0 * ^€=0 
a Vo = %=0 + % 

where 


H 

^0 


= measured dynamic amplitude 
= incremental error 

= measured dynamic amplitude with error 
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and 


a-i - measured phase angle 

e= 0 

a-i - error In a-, 

e ± 

a. , = measured phase angle with error 

i e^0 


e„ = 0.01, 0.02, 0.05, 0.04, 0.05 

9 



Equation (42) is solved for the aerodynamic derivatives based on 

the response measurements with errors. These results are presented in 

Table V. Amplitude errors of the form are normally associated with 

instrument calibration errors, numerical results show that does 

not affect C 7 and appears to have a one to one effect on C, . 

l p t 5 

is associated with errors due to the sin-cosine resolver. Phase 
angles determined with this equipment were repeated several times at 
each frequency, and the repeatability was normally within 1° to 5° • 
Numerical results show that phase errors affect C 7 slightly and C 7 

h p 

about 2 per cent per degree of error depending on the dynamic pressure. 
In general, this simplified analysis shows that realistic values of 
response errors (e^ < 0.05, cu^ < 5°) predict results within about 
5 per cent accuracy. 

4.3.5 Comparison of Experimental Results 

A comparison of the aileron derivative C 7 , measured both 

1 5 

statically and dynamically, is presented in Figure 12. Results are 


presented in terms of 



versus model dynamic pressures at Mach 


TABLE V.- ROLL RESPONSE ERROR ANALYSIS 


0 C Z p / rad 


-0.310 

-0.310 

-0.310 

-0.310 

-0.310 

-0.510 


-0.29k 

-0.29k 

-0.29k 

-0.29k 

-0.294 

-0.294 


0 .0123 0 


a i e > 

de| 

C 7 /rad 
P 

Ci /rad 
L b 

0 

-0.3^9 

0 .0166 

1.0 

-0.356 

0.0168 

2.0 



3.0 

-O.369 

0.0171 

4.0 

-0.375 

0.0173 

5.0 

-0.381 

0.017^ 

0 

-0.310 

0.0123 

1.0 

-0.313 

0.0123 

2.0 

-0.316 

0.0124 

3.0 

-0.318 

0.0124 

4.0 

-0.319 

0.0124 

5.0 

-0.321 

0 .0124 I 


0.0107 

3.0 

0 .0108 

4.0 


-0.294 0.0104 I 0.675 I 150 

-0.300 0.0105 
-0.306 0.0106 
-0.311 0.0107 

-O.317 0.0108 

-0.321 0.0109 
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numbers 0.375 and 0.75. Note the close correlation between static and 
dynamic results at all but the higher dynamic pressures. These results 
tend to confirm the assumption made earlier: that as the ailerons 

become less effective, the dynamic motion can no longer be approximated 
by a single-degree-of -freedom solution. Since C ^ cannot be measured 
statically, experimental data are not available for a comparison of the 
damping -in -roll derivative. 

The reversal boundary was also determined statically at Mach 
numbers 0.675 and 0.75 using the two-cable mount. These results are 
presented in Figure 9 to give comparison of static testing procedures 
using both the sting-pylon-spring and two-cable mounts. Note the 
extremely close comparison of results between the two mounts for deter- 
mining reversal dynamic pressure. 


CHAPTER V 


SUMMABY AND CONCLUSIONS 


A new wind-tunnel technique for measuring various aerodynamic 
derivatives of an aeroelastic model has been presented. The technique 
applies free-flight procedures to a model flown in the wind tunnel on 
the two-cable-mount system. The complete equations of motion have been 
derived in terms of model properties, mount geometry, and aerodynamic 
derivatives . 


In the case of the longitudinal equations of motion, it is theo- 
retically possible to uniquely determine each of the aerodynamic deriva- 
tives by measuring the model response to a steady-state sinusoidal 
oscillation of the horizontal tail. In the determination of the deriva- 
tives from model test data a least-squares procedure is used to solve 
the set of redundant equations generated. A numerical example has shown 
that the derivatives can be determined uniquely only if exe~t response 


data are analyzed. The derivatives Cj^, C^, C m ^ , and (C^ + ) 

can be predicted with significant accuracy using this technique. An 


alternate static method to measure has been presented. The assump- 

CmR c 

tion that Cj^ = — provides an estimate for C^. A comparison 
between flight and wind-tunnel equations shows that, due to the added 


mount system restraints, the equations can be solved for each of the 


derivatives uniquely. However, introduction of an error into the model 


response causes the solution to become ill-conditioned resulting in 
equations quite similar to those used for determining the aerodynamic 
derivatives in free flight. 
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In the lateral equations of motion a basic free-flight assumption 
of single-degree-of -freedom response in roll allows the experimental 
verification of the dynamic approach to derivative measurements . The 
model in the wind tunnel can be forced to behave in a similar fashion 
by the proper selection of mount system parameters. The derivatives 
Cjg and can be determined by measuring the dynamic model response 

to a steady-state sinusoidal oscillation of the ailerons. Experimental 
results obtained on a l/l9-size aeroelastically scaled model, tested 
both statically and dynamically in the wind tunnel to determine C^, 
verifies the application of this new testing procedure. 

It has been shown that by means of a rather simple two -cable -mount 
system the dynamic characteristics of an aircraft can be closely simu- 
lated, allowing the use of free-flight techniques to estimate the aero- 
dynamic derivatives of an aircraft in the early design stages. Application 


of the analysis presented in this thesis should assist in developing 


testing techniques required to satisfy specific research programs. 
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APPENDIX A 

MOUNT RESTRAINT INFLUENCE COEFFICIENTS 
Vertical Translation 

As shown in Figure 3 the cable configuration analyzed has a 
vertical forward cable and horizontal rear cable. The linearized mount 
stiffness due to small perturbations about a trim point will be deter- 
mined. A diagram of the forward pulley configuration is shown in 
sketch (1A) . 



Sketch (1A) 
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The Z force at each pulley due to a z displacement is 


Z 1 = -T F sin(P F + Ap x ) + T f 
Z 2 = T f sin(P F + Ap 2 ) - Tp 


(A-l) 


Since all displacements are small perturbations about the trim point 


Z 1 + Z 2 = T f (A 3 2 - A3 1 )cos pp (A-2) 

AP-^ ai;d Apg are composed of angular changes due to displacements in 
the x and z directions . 


x displacement: 



xi << Lp 

From the geometry 


x 1 sin P F sin(p F + APp) 
sin(p F 4 Ap x - p F ) 

x 1 sin Pp sin(3p + APj) 
sin(AP 1 ) 


For -'mail perturbations 


sin AP^ = AP X 
sin(P F + Ap^ = -^ 

Jjp 


hence, 


APi = sin pp 
T 1 


z displacement: 



From the geometry 


Z 1 < < z o 


x sin p F sin(p F + AP 2 ) 
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For small perturbations 


sin(3 F + AP X ) = 


z o + 2 1 

Lp 


hence 


Ap l ~ cos % 


(A-4) 


Therefore, the total A3^ due to an x and z displacement is 


- 


sin pp + z^ cos pp 


(A-5) 


In a similar manner it can be shown that 


AP 2 = 


x 2 sin pp - z 2 cos Pp 


(A-6) 


Pulley movements due to z and 0 displacement (sketch (1A)) are as 
follows : 


Xj = -h9 


x 2 = h0 


z-^ = z - e0 


Zo = z - e0 


(A-7) 


Substituting equations (A-5)# (A-6), and (A-7) into (A-2) gives 


Z-^ + Z2 “ 


2Tt? — 

— cos Pp £he sin Pp + z cos Pp 


- e0 cos 


h3 


(*- 8 ) 
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A diagram of the rear pulleys is given in sketch (2A) 



Sketch (2A) 

The vertical force at each of the rear pulleys due to a z displacement 
is 




(A-9) 



2T r 

z 3 + z b = ~ L^~^ Z + a9 ^ 


(A-ll ) 
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The sum of the vertical forces at each of the pulleys ((A-ll) and (A-8)) 
due to a z and 0 displacement is 


Z c — Z 1 + Z 2 Z^ + Z|^ 


Zc = -2zU- cos^pp. + 




h cos P F sin 0p 


i v 2 

e cos dp - a 


(A-12) 


Since 0 = 0 + 0 t , we can express equation (A-12) in stiffness-influence 


coefficient form as 


Z C ~ Z C D " K ZZ z " W 


(A-13) 


Where 


K zz = 2 [if cos % 


' + *9 


(A-lfc) 


Td Tp Tp o "1 

K^ e = 2 a — - h cos 0 F sin 3 F - e -*■ cos^(3p 

L Lr 4 1 rF J 


(A-15) 


Zc 0 = -K ez 0 t 


(A-16) 


Pitch 


Tlie fortes producing moments about the center of gravity from the 
forward cable are (sketch (1A)): 
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z x = “Tp sin(Pp + Apj) 
Zg = Tp sin(j3p ^ A3 2 ) 
X x = Tp cos ( 3p + A3^ ) 

= Tp cos(pp + A 32 ) 


(A-17) 


For a positive z, 0 displacement, the moment produced by the forward 
cable is 

Mq = -Z x (e - he) - Z 2 <e + he) - X^h + ee) + X 2 (h - e?) 

F 

(A-18) 

Substituting equations (A-5), (A-6), (A-7)> and (A-17) into equa- 
tion (A-l8) and simplifying gives 


Mq = 2Tp9 -h sin Pp - e cos Pp - 


(e cos pp + h sin Pp)‘ 


2 Tpz — p 

+ ~ |e cos^Pp + h sin Pp cos PpJ 


The forces producing moments about the center of gravity from the rear 
cable sure (sketch (2A)): 
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(A-20) 


-T r cos 3 r 


-T r cos 0 r 


The pulley displacements due to a positive z, 9 displacement are 


= z 4 = z + a9 


x, = x,. = a 


3 - 


The moment about the center of gravity produced by the rear cable is 


Mq = Z^a + Z^a + XjaG + X^a0 


Substituting (A-20), (A-21), into (A-22) and simplifying, results in 


“Cjj “ 29 “ “j““ " aT R cos PR 


“c - % + \ 


The total moment 


101 


Hence, 


— a 2 Tn f” 1 

Mq = 29 - -£-£• - aT R cos Pp - Tp h sin 3p + e cos Pp + ^(e cos Pp 

+ h sin Pp)^| + 2z|j + --(e cos 2 Pp + h sin Pp cos pp)j 


As before, 9 = 6 + 0^. and -we can express equation (A-24) in stiffness- 


influence coefficient form as 


Where 


M c - " K dz z " ^e 9 


Mcq “ _K ee 9 t 


2^|p 2T 

Kez = Ijf ’ COB \ + h sin h cos M 


2a 2 T R f 

Kge = — + 2aT R cos P R + 2T p h sin Pp + e cos Pp 


+ i-(e cos Pp + h sin Pp)^j 


Note that for small perturbation theory K^q = Kq£. Equations (A-13) 
and (A-25) completely define the mount system restraints in terms of 
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the form 



h: = _K ^Y y " ~ 

N c = -K^yy - - K^'lf 

Y C = " K YY y “ K Y0^ “ 
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Where 


K YY * + ^ COB % 



0 


*» * — - V 006 M d 8ln 


Pr 


% - k Y5( “ 0 

*# - ai p(^ * 8ln hj + 2 ®r(4 + 8i " %) 
v*-° 

^ “ ity “ 0 


Kyy = Kyy 


Kyy * 2eT] 


^ + C08 %) + 2t r(® 

+ ^(a cos Pp d sin 3 r) 2 J 


a cos Pp + d sin Pp 


APPENDIX B 


LEAST-SQUARES SOLUTION OF REDUNDANT LINEAR EQUATIONS 
General Solution 

Given at set of N equations in M unknowns 

Y 1 = A ll x l + a 12^2 + a 13 x 3 + • • • Alm^n 

Y 2 = A 21 X 1 + A 22 X 2 + ^ 3 X 5 + . . . A 2 m X m 

% = A N1 X 1 + A N2 x 2 + a N 3 x 3 ■*■••• a NM% n > M 

Equations (B-l) can be written in the form 

M 

h ■ ) Vi. 1 - 1 > 2 ’ 3 ' • • • " < B - 2 > 

m=l 


We wish to solve for (X^, Xg, . . . X m ) given (Y^, i = 1,2, . . . Nj 
A im» m = • • • M) . Choose (X^, Xq, . . . 5^) so that the sum of 

the squares cf the deviations is as small as possible. Where 

M 

v i - £ <*iA - Y i> 

in=l 

is the deviation. Therefore, 

N 

li- ^ (Vt ) 2 
i=l 



10b 
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is to be minimized. Hence, 



where 


k = 1,2,5, . . . M 




m = k 
m / k 



Equation (B-5) results in M equations which are solved simultaneously 
for X^, X 2 , Xj, . . • Xjjj. 
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Complex Solution 

If the Y's and A's are complex, define 


N 


N 


*- > N - > 


1=1 


i=l 


M 


^ 

m=l 


where 


Yi = Yi + lY t 
x •‘■real imag . 


^im A i® re al + ^^“imag. 


also, the complex conjugates of these quantities are 


Y* = Y, - 1Y, 
x x real imag . 


A lm A im rea i ” ^imimag. 


(B-4) 


Therefore, equation (B-4) is written 

N / M 


F M 


E 


- ) ) A im*m - Yi f ) ^ 


i=l\m=l 


\m=l 


which is to be minimized. Therefore, 

Be 


Sx k 


= 0 


k — 1,2,3, • • 
m = k 


. M 


m 4 k 
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Hence, 


N 


V 

1=1 


M 


M 


), A i^j A ^ + ( ^ 


N 


- ) <n A it * n A i k ) 


i=l 


(B-5) 


Equation (B-5) results in M real equations which are solved for 

^L> ^2> Xj, . . . . 


Example: M = 5 

k = 1 

N 


C A il X l + A i2 X 2 + A i3 X 3 )A*i + (AJ 1 X 1 + A^Xg + A^X^aJ 


N 


" ) ( Y i A !l + X iA?l) 


i=l 


2X, 


a, a, n 

L M 8 * *2 I (A Ii A 12 - A 11 A 12> + *J ) ^ A ll A i3 + A 11 A ;,) 

1=1 1=1 U 1 


N 


• ). (I i A u + Vh> 


i=l 


Similarly for 
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APPENDIX C 


COMPUTER PROGRAMS 
Model Response 

Assuming sinusoidal motion in the form z(t) = , 

9(t) = equations (17a) and (17b) can be written. 

Vertical translation: 


z Q e 


L$l 


{-a?m + K zz ) + ioj^Cr^ + . D ) 


+ 0 o e 1 ^ 2 qSC Iti = -qSC^5 0 


Pitch: 


i01 qSc 2 o£ n ^ qSc 

_2 


z o e rx | — ' ~ C m - ico^ c 

,2U 2 “ U 


+ - qScC^ + 


qSc 2 


- *“ + Vj ‘ '’ C? V° 


The equations above can be written in the form 

;1 ^[ A 1 + iB l] + = F c 


z e 
o 


(c-1) 


+ IBj] + + iB 4 ] = 


M„ 


(C-2) 
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Where 


A^ = -a?m + 1(22 

A qSc 2 o£ 

" 5= 2U 2 ^ 

B 1 = ® + c D ) 

b 3 = -irSx 

Ag = QSCj^ 

A k = -1^ _ qScC % + ^ 

F o= ^ SC Lb 5 o 

_2 

+ V 

Applying Kramer's rule results in 

»c • 


z o e 


% _ 


Fo A 2 

A4 + iBij. 
Dp + 1% 



A 1 + % F 0 

A 3 +1B ? -o 
Br + iD i 


Where 


Dr - A^Ajj. - AgAj - B^Bl). 


D i = AiB,. + - AgBj 


Ill 


Each of the terms z^ 1 ^ 1 and 9 c> e i ^ 2 can then be divided into the 
real and imaginary parts z Q , e o , , and • A i> A 2» A 3> A4, 

Bj_, Bj, and are defined in the program listing. Once the values 

of the aerodynamic derivatives, model physical properties, and test 
conditions are selected equations (C-l) and (C-2) are solved simultaneously 
to determine the dynamic response z Q , 0 Q , 0^, and as a function 

of cu. A sample program listing is presented. 


Definition of terms : 

M 

- mass of model, slug6 

KZZ 

- vertical spring restraint, lb/ft 

Q 

- dynamic pressure, psf 

s 

O 

- wing area, ft 

u 

- wind-tunnel flow velocity, ft/sec 

c 

- mean-aerodynamic chord, ft 

IY 

- pitch inertia, slug-ft^ 

CIA 

- lift curve slope Cj^/rad 

DC 

- drag coefficient Cjj 

CMADOT 

' c m A / rad 

CMA 

- C^/rad 

CMQ 

- Cn^/rad 

FO 

- _ qSCj^6oj lb 

MO 

- qScC^bo, ft -lb 

N 

- number of frequencies at which response is calculated 

WW(I) 

- values of cd to be used to calculate response 
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Where 

B(l,l) 


B(l,2) 


B(l,3) 


B(l» 


A1 - A 1 

A2 - Ag, etc. 


Error Analysis 

Equations (20) and (21) are programed in the following manner: 


— 

B(l,l) 

B(l,2) 

B(l,3) 



B(l,4) 

B(2,l) 

B(2,2) 

B(2,3) 

c D 

= 

B(2,4) 

B(3,l) 

B(3,2) 

B(3,3) 

% 


B(3A) 


N 


l 

J=1 


N 

. .2 n 

|Aji| 2 = (qS) 2 £ 

O /^r. 1 2OB„0„ 

60 + ( u J + D 8ln( ^ ‘ ^ 

1-1 

- _ 


OPCOi 


H, o N r 

= \ (A * A + A A* ) - S - (qS) V {a * oY 

L y V J 2 A J1 j2 ' u ) I u 

J=1 J=1 


J=1 

+ 0 » o e o sin(0 2 - 0i) 


O>=0D < 


N 

l 

J=1 


^ A ji A j3 + a ji a j3^ ' 


2q2 S 2 6 o ^ fe o cos ^ ^ sin 

J^l L J 


aj»G0 4 


N, N 

I (b J A Jl + b J A Jl } “ 2<JS Z ^ " K ZZ )z o 6 o COs(< *2 " ^lO 
J=1 J=1 


o*=a>4 


115 


Where 

A(l,l) = 


A(l,2) = 


A(l,3) = 


A(1,0 = 


A(l,5) = 


A(l,l) 

A(l,2) 

A(l,3) 

A(1A) 



A(l,5) 

A(2,l) 

A(2,2) 

A(2,3) 

A(2,4> 

Si 


A(2,5) 

A(3,l) 

A(3,2) 

A(3,3) 

A(3,4) 

Si 


a(3,5) 

a(M) 

A(4,2) 

A(4,3) 

A(4,4) 

f“6_ 


A(4,5) 
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The program generates the coefficients A(i, j ) , B( i,J ) based on 
the test conditions, model properties, and model dynamic response. The 
effect of an error in response is introduced by considering 


Vo 

Vo 

e^O 

Wo 


= z (1 + Ac z ) 
e=0 

= 0 (1 + A€ 0 ) 

°e=0 w 

' ^.0 + 

‘ %,o * 


Next, the program solves equations (C-3) and (C-4) separately for the 
value of each aerodynamic derivative with response errors included. 
(The equations are solved by a library subroutine called MATRIX.) 

A sample program listing is presented. 

Definition of terms: 

U - wind-tunnel flow velocity, fps 

Q - dynamic pressure, psf 

p 

S - wing area, ft 

M - model mass, slugs 

C - mean aerodynamic chord, ft 

IY - pitch inertia, slug -ft ^ 

KZZ - vertical mount restraint, lb/ft 

KTT - pitch mount restraint Kq 0 , ft-lb/rad 

DELT - tail angle 6 Q , rad 

ZE - Ae^ 
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TE 

" A€ e 



PIE 

- A< %. 



P2E 




K 

- program code 0 

new test case 



1 

new error function 



2 

stop 

N 

- number 

of response data to be analyzed 

W(I) 

- values 

of tail frequency 

zo(i) 

- values 

of response z Q 


THEEAO(I) 

- values 

of response 0 O 


PHIl(l) 

- values 

of phase angle 

h 

PHI2(I) 

- values 

of phase angle 

$2 

MATRIX 

- solves 

equations (C-5) 

and (C-4) 
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MODEL RESPONCE PROGRAM 


DIMENSION IN (100) 

REAL KZZ,KTT,M,MO, I Y 

1 F0RMAT(10F8.0) 

READ ( 1 ) M/ KZZ/Q, S, U, C, I Y,KTT, CLA, CD/ CMADOT, CMA, CMQ 
PRINT 2 

2 F0RMAK25M I ABEL ‘AIL STOP 340//) 

REA0(5/1) FO/MO 

PRINT 3/FO/MO 

3 FORMAT (/ 511 FO -F10.4/5IJ MO -F10.4//) 

PRINT 4 

4 FORMAT (16H OMEOA - RAD/SEC/ 8X/ 7IIZ0 - FT, 10X, 12HTHETA0 - RAD, 
16X/10HPHI1 - DEG/6X/10HPHI2 - DEG//) 

5 FORMAT (14) 

6 FORMATUOF8.2) 

READ(5, 5)N 

READ(5,6)(UW(I ),I-1,N) 

1-0 

7 1-1*1 
W-WW(I ) 

A1-KZZ-M*W**2 
Bl»(CLA*CD)*W*n < :/lJ 
A2*Q*S*C LA 

A3*Q*S*C**2*W**2*C‘1ADOT/2. 0/U**2 

1(3— W*Q*S*C*CMA/U 

\4— IY*W**2*KTT-Q*S*C*CI1A 

B4— (CMAOOT*CMQ)*W*Q*S*C**2/2 . /U 

DR-A1*A4-A2*A3-B1*B4 

DI«A1*B4+A4*B1-A2*B3 

ZO-SQRT((DR*(FO*A4-MO*A2)*DI*FO*B4)**2*(OR*FO*B4-DI*(FO*A4-MO*A2)) 
1**2 )/ »DR**2*DI **2 ) 

PHI 1-57. 295780*ATAN2(DR*FO*B4-DI*(FO*A4-MO*A2),Dk*<FO*A4-MO*A2)*DI 
1*F0*B4) 

THETAO-SQRT((DR*(MO*A1-FO*A3)*DI*(MO*B1-FO*B3))«*2*'.DR*(MO*B1-FO* 
1B3)-DI*(M0*A1-F0*A3))**2)/(DR* 2*DI**2) 
PIII2«57.295780*ATAH2(DR*(M0*B1’F0*B3)-DI*(M0*A1-F0*A3)/DR*(M0*A1- 
1F0*A3)*DI *(M0*B1-F0*B3) ) 

PRINT 8/W/Z0/TIIETA0, PHI 1, PHI 2 

8 F0RMAT(4XF6.2,8X4E18.8) 

PUNCH 0, W, ZO, THF.TAO, PH 1 1, PH I 2 

9 FORMAT(5E15.8) 

IF(I.LT.N) GO TO 7 
STOP 

END 
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ERROR ANALYSIS PROGRAM 


DIMENSION W(50)/Z0(50)/THETAO(50)/ PHI 1(50), PHI 2(50) 

DIMENSION A(4,5),B(3,4) 

REAL KZZ,KTT,M,IY 
1 FORMAT(5F15.8) 

3 F0RMAT(30X,42I! LONG. LEAST SQUARES SOLUTION - TAIL INPUT////) 

4 FORMAT (511 H-F10.4,4X, 511 Q-F10. 4, 4X, 511 S«F10.4,4X, 5H M-F10. 

14/ 4X/ 511 C-F10.4//5H IY-F10.4,4X,5H KZZ-F10. 4, 4X, 5H KTT-F10.4,4X 

2,5HDELT«F10.6,4X,5H ZF.-F10.4//5H TE-F10. 4,4X, 5H P1E-F10.4/4X, 511 

3P2E-F10.4////) 

5 FORMAT (14) 

6 FORMAT(5E15.8) 

PRINT 3 

20 READ(5,5)K 

IF(K.EQ.O) GO TO 14 
IF(K.EQ.l) GO TO 16 
IKK.EQ.2) GO TO 1500 

14 READ (5/ 5 )N 
DO 15 l-l/N 

15 <<EAD(5,G) W(l )/ZO(l )/THETAO(l )/PMIl(l )/PHI2(l ) 

16 READ(5/1) U,Q,S,M,C, IY,KZZ,KTT,0ELT,ZE,TE,P1E,P2E 
PRINT 4/U,Q/S/M,C,IY,KZZ,KTT,DELT,ZE/TE/PlE/P2E 
N-l 

d-30 

11 DO 101 11-1/4 
DO 101 Jl-1/5 

101 A( I 1/J1)«0. 

DO 102 11-1/3 
DO 102 Jl-1,4 

102 B( 1 1/ Jl)-0. 

PRINT MM 

9 FORMAT (3H N-l 2, 6X, 2HJ-I 2/// ) 

DO 100 l-N/J 

Z0(l )-(1.0*ZE)*ZO(l ) 

THETAOO )-(1.0+TE)*THETAO( I ) 

PHI 1(1 )-PlE*PH'l(l) 

PH 12(1 )-P2E+PHI 2 ( I ) 

A(1/1)-A(1/1)*(W(I )*THETA0( I ) )**2 

A(1,2)-A(1,2)*(THETA0(I)**2*W(I )*THETA0(l )*Z0(l )/U*SI N ( (PH I 2 ( I )-PH 
111(1 ) )/57. 2958 ) )*W( I )**2 

A(l/3)-A(l,3)*W(l )**2*THETA0( I )*Z0( I )*COS( (PH 11(1 )- PH 12(1 ))/57.295 
18) 

A(1,4)-A(1,4)*W( I )*THETA0( I )*S I N(PHI 2( I )/57.2958) 

A(2, 2)-A(2, 2)+(THETA0( I )**2*(W( I )*Z0( I )/lJ)**2*2. 0*THETA0 (I ) *Z0 ( I ) * 
1W(I )/U*SIN{(PIII2(l)-PHIl(l ))/57.2958))*W(l )**2 
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A(2,4) a A'2,4)*W(l )*THETAO( I )*S IN(PHI 2( I )/57. 2958)+W( I )**2*Z0( I )/U* 
1C0S(PIII1(I )/57. 2958) 

A(2 / 5)“A(2 / 5)*(IY*W(I )**2-KTT)*W( I )**2*ZQ( I )*THETAO( I )/U*COS( (PHI 2 
l(l)-PHIl(l))/57.2958) 

A(3,3)>A(3,3)*THETA0(l)**2*(W(n*Z0(l >/U)**2+2.0*THETA0(l >*ZO(l )* 
1W(I)/U*SIN((PHI2(I )-PHIl(l ))/57.2958) 

A(3/4)“A(3 # 4)*THETA0( I ) *COS(PHI 2(1 )/57.2958)-W(l )*Z0< I )/ll*SIN(PHIl 
1(1 )/57. 2958) 

A(3,5)«*A(3,S)+(IY*U(I )**2-KTT)*(THETA0( I )**2*U(I )*7.0(l )*THETAO(l )/ 
1U*SIN( (PHI 2(1 )-PHI 1( I ))/57.2958)) 

A(4,4)-A(4,4)*(Q*S*C*DELT)**2 

A(4,5)-A(4,5)*(IY*W(I )**2-KTT)*TIIETA0(l )*C0S(PHI2(I )/57.2953) 
B(1 / 1)*B(1,1)+TIIETA0( I )**2+(U(l )*Z0( I )/U)**2*2. 0*W( I ) *THETAO( I )* 
1Z0( I )/U*SIN((PHI2(l)-Pllll(l))/57.2958) 

B(1,2)-B(1,2>*W(I )**2*Z0(I)**2/U*W(I )*THETAO(l )*Z0( I )*SIN( (PHI 2( I ) 
1-PIIIKI ))/57.2958) 

B(1,3)«B(1>3)+THETAO(I)*COS(PHI2(I )/57.2958)-W( I )*7.0( I )/U*5 1 M(PHI 1 
1(1 )/57. 2958 ) 

B(l,4)«B(l,4)+(tl*W(l )**2-KZZ)*7.0(l )*THETAO(l ) *COS( (PHI 2( I ) - PH 1 1 ( I ) 
D/57.2958) 

B(2 / 2)“B(2/2)+(W(l )*ZO(l ))**2 
B(2,3)»B(2,3)+W(I)*Z0(I )*SIN(PHI1(I )/57.2958) 
B(3,3) a B(3,3)*(Q*S*DELT)**2 

B(3,4)«B(3,4)*(H*W(I )**2-K7.Z)*Z0(l )*C0S(PHI1(I )/57.2958) 

Z0( I )*»Z0( I )/ (1.0+ZE) 

THETAOd )»THETAO(l )/(1.0+TE) 

PHI 1(1 )-PHIl(l J-P1E 
100 PH 12(1 )=PHI 2( I )-P2E 

ACl,l) a A(l,l)*(Q*S*C**2/2.0/U)**2 

A(l,2) a A(l,2)*(Q*S*C**2/2.0/U)**2 

A(l,3) a A(l,3)*(Q**2*S**2*C**3/2.0/U**2) 

A(l,4) a A(l,4)*(-Q**2*S**2*C**3*DELT/2.0/IJ) 

A(2/l) a A(l/2) 

A(2,2) a A(2, 2 )*(Q*5*C**2/2 . 0/U)**2 
A(2/4) a A(2/4)*(-Q**2*S**2*C**3*DELT/2.0/IJ) 

A(2/ 5 ) a A(2, 5)*Q*S*C**2/2. 0/U 
A(3 / l) a A(l / 3) 

A(3,3) a A(3,3)*(Q*S*C)**2 

A(3/4)“A(3/4)*(Q*S*C)**2*DELT 

A(3,5)»A(3,5)*(-Q*S*C) 

A(4,l) a A(l,4) 

A(4,2) a A(2,4) 

A(4,3) a A(3,4) 

A(4,5) a A(4,5)*(-Q*S*OOELT) 

B(l,l) a B(l,l)*(Q*S)**2 
B(l,2) a B(l,2)*(Q*S)**2/U 
B(l,3) a B(l, 3)*(tQ*S)**2*DELT 
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B(1,4)-B(1,4)*Q*S 

B<2,l)-B(l,2) 

B(2,2)»B(2,2)*(Q*S/U)**2 

B(2,3)-B(2,3)*(-Q*S)**2*DELT/U 

B(3,l)-B(l,3) 

B(3,2)«B(2,3) 

B(3,4)«B(3,4)*(Q*S*DELT) 

PRINT 400 

400 FORMAK//9H A MATRIX//) 

PRINT 8,((A(ll,Jl),Jl a l,5)s 11*1/4) 

8 FORMAT(5Elfi.8) 

PRINT 401 

401 FORMAT(//9H B MATRIX//) 

PRINT 10,((B(I1,J1),J1«1,4),I1*»1,3) 

10 FORMAT (4E1B. 8) 

CALL MATR I X ( 10, 4, 5 , 0, A, 4, DETERM ) 

CALL MATRIX(10,3,4,0,B,3,DETERM) 

PRINT 7,A(1,5),A<2,5),A(3,5),A<4,5),B<1,4),B(2,4),B(3,4) 

7 FORMAT ( 5H CMQ-F10.4//8H CMADOT-F10.4//5H CMA-F10.4//7H CMDEL-F10.4 
1//5H CLA-F10.4//4H CO-F10.4//7H CLDEL-F10.4/// ) 

I F ( J.GE. 30) GO TO 20 
1500 STOP 
END 



